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ABSTRACT 

This work further develops the properties of fractional differential forms. In particular, 
finite dimensional subspaces of fractional form spaces are considered. An inner product, 
Hodge dual, and covariant derivative are defined. Coordinate transformation rules for 
integral order forms are also computed. Matrix order fractional calculus is used to define 
matrix order forms. This is achieved by combining matrix order derivatives with exterior 
derivatives. Coordinate transformation rules and covariant derivative for matrix order forms 
are also produced. The Poincare lemma is shown to be true for exterior fractional 
differintegrals of all orders excluding those whose orders are non-diagonalizable matrices. 
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1. Introduction. 

The use of differential forms in physics, differential geometry, and applied mathematics 
is well known and wide spread. There is also a growing use of differential forms in 
mathematical finance [23], [3]. Many partial differential equations are expressible in terms 
of differential forms. This notation allows for insights and a geometric understanding of 
the quantities involved. Recently Maxwell's equations have been generalized using 
fractional derivatives, in part, to better understand multipole moments [5]-[8]. Kobelev has 
also generalized Maxwell's and Einstein's equations with fractional derivatives for study on 
multifractal sets [10]-[12]. Many of the differential equations of finance have also been 
fractionalized to better predict the pricing of options and equities in financial markets [1], 
[22], [24], [25]. Fractional calculus is also being applied to statistics [18]-[20]. Hence, 
there is a need to further develop the notion of a fractional differential form to better 
understand fractional differential equations, fractional line elements, and fractal geometry. 

In an earlier work [4] (hereafter referred to as FDF I) the study of fractional differential 
forms was initiated. The purpose of that study was to combine Riemann-Liouville fractional 
calculus with exterior calculus on an n-dimensional Euclidean space. FDF I focused 
primarily on differential order fractional forms. In this paper the properties of differential 
order fractional forms are extended. Integral and matrix order forms are also introduced 
and studied. Following Oldham and Spanier [17] the term fractional differintegral form 
shall be used when the result applies to both differential and integral order forms. 

In section 2 some algebraic properties are worked out for various finite dimensional 
subspaces of the F(v,m,n), vector spaces defined in FDF I. In particular, a Euclidean inner 
product and Hodge dual are constructed for these finite dimensional subspaces and then for 
the infinite dimensional F(v,m,n) . In section 3 some differential properties are worked out 
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for the finite dimensional subspaces of F(v,m,n) . The fractional Poincare lemma is 
considered for all orders of the fractional exterior differintegral, including orders that 
represent fractional integration. Coordinate transformation rules for integral order forms 
are constructed. An example of a change of coordinates for integral order forms is also 
provided. In section 4, matrix order differintegrals are used to define matrix order forms 
(see the appendix for the definition of matrix order derivatives and integrals). Coordinate 
transformation rules and a covariant derivative are also constructed. The fractional Poincare 
lemma is extended to matrix order forms and found to be true provided that the matrix 
representing the order is diagonalizable. 

As in FDF I the coordinate index will be a sub-script rather than the traditional super- 
script. The summation convention will not be used so as to maintain some clarity with the 
myriad of indices to be encountered. Throughout this paper shall represent Cartesian 
coordinates and, {yj and {z,.} shall represent arbitrary curvilinear coordinates on an n- 
dimensional Euclidean space, E" . a i shall denote the initial point for evaluation of the 

fractional differintegrals in Cartesian coordinates and, a i and a t in the curvilinear 
coordinates. 

Before beginning it is necessary to correct an error and two typographical errors from 
FDF I that may cause confusion. In FDF I equation (25) is in error and equations (26) and 
(27) have miss typed indices. The equations should read, respectively, 

{dx^ Adx? | i,j e {1,K ,n}y^ 2 > 0, 3 fl, + fl 2 = v] 

p=£ £ j v (Av(v,v-v,.)^r A^p) jv, 

i=l 7=1 
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P = ZZZ\ V ^""'(A^CVpV-v^v-v.-v,)^' *dx]> Adxr^dVjdv, 



(3) 



A more compact notation for fractional differintegrals is also adopted. The notation for 
fractional integrals and derivatives is given below. 



(4) 



1 \ J a . 



A-n+1 



r(n-A) Jfl <(*-§) 



Re (A) > 



n > Re(A) (n is whole) 



(5) 



The parameter X is the order of the integral or derivative and is allowed to be complex. In 
the appendix it is shown that % may assume matrix order values as well. Equation (4) is a 
fractional integral and equation (5) is a fractional derivative. Taken together, the operation is 



referred to as differintegration (see page 61 of [17]). 



Equations (6) - (13) are the basic identities from fractional calculus. In the following, m 
is a whole number and p and q are complex numbers whose real part is greater than zero. 



dx' 



(6) 



a D^D-f(x) = f(x) 



(7) 
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. P7,a D'fW * f(x) 



(8) 



D-f(*) = i( D»fW (9) 



D p D q f(x)= D p+9 f(x)-E D ,_y 'f(jc)| ^' aJ 

fl. x, a, x, V / a, x. A V-^v a, x f A V-^/ T -,/ 1 .\ 

r ('-/ ; -/) 



(10) 



„ D- fl| D'f (*)=., DxT'f W -t a, V%%x)\ ^'° ,)P [ (ID 

; =i k=fl ' r(l+p-;) 



a D;"( fl D-(fW)) = a D^>(f(*)) (12) 



s=0 



V s J 



(13) 



In equations (10) and (1 1) k is the first whole number > Re(q) . The above formulae and 
definitions can be found in [16], [17], [21]. 

2. Algebraic Properties. 

In this section the restriction that the order of the fractional differintegrals be real is 
made. Many of the results will clearly be valid for complex order differintegrals as well. 
To begin, define a vector space G(v,n) that is a finite dimensional subspace of F(v,m,n) 
(this is defined in FDF I) such that, 
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m 

G(v„v 2 , K ,v m ,n) = nG(v,,n), 



(14) 



and, 



F(v,m,n)= U G(v,,v 2 , K ,v m ,n). 



(15) 



If V > then v ; > for all i is imposed. If V < then v i < for all i is imposed. 



n 

cc g G(Vj,v 2 , K ,v m ,n)^a= J^a. L (m dx^ aL acU 



(16) 



Some of the terms in the sum may be zero if any of the orders of the coordinate differentials 
are the same. For example, dx 1 / 2 Adx]' 2 = 0. The dimensions of these vector spaces are 



dim(G(v p n)) =n, 



(17) 



dim(G(v 1 ,v 1 ,n)) = 



i \ 
n 

V 2 / 



(18) 



dim(G(v 1 ,v 1 , K ,v p n)) = 



i \ 
n 



m<n. 



(19) 



In equation (19) there are m coordinate differentials of order v r 



dim(G(v 1 ,v 2 ,n)) = n v 2 



(20) 
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dim(G(v„V 2 , K ,v m ,n)) = n 



(21) 



dim(G(v 1 ,K v 1; v 2 ,K ,v 2 ,K ,v m ,K v m ,n)) = 



n n n 

L 

\Pl)\Pl) \PmJ 



(22) 



In equations (21) and (22) the v, are distinct for i e {1,2,K ,m} . In equation (22) there are 
p l coordinate differentials of order v v p 2 coordinate differentials of order v 2 , etc., i.e. 



m P i 

G(v„K v„v 2 ,K ,v 2 ,K ,v m ,K v m ,n) = n ri G ( v /' n 



7=1 



(23) 



An inner product on G(v p n) can be defined in the usual way (see [2], [9], and [14]) and 
then built up to an inner product on the higher dimensional subspaces. Let a, j5 e G(v 1? n) 

n n 

with a = ^a i dxj l and /? = 52/? 1 d* ; I Vl • The Euclidean inner product on E" for G(v l5 n) is 

i=l i=l 



(24) 



In curvilinear coordinates equation (24) becomes 



(25) 



',7=1 



The metric ^ y (y,Vj) is defined in equation (67) of FDF I. The cbc ( v ' form an orthonormal 
basis on G(v l5 n). 
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(26) 



If it is assumed that the v,. are distinct for i e {1, K , m} , then the inner product on 

n 

G(v l5 v 2 , K ,v m ,n) is given as follows. Let a = £a,. L ; cbc, Vl aL AdxJ'" and 



= EA l im dx,-' AL Adx^then, 



'l L 'm '] 

I.L i" =1 



(«,£)= • (27) 



/jL l m ~ l\ 

L L i" = 1 



In curvilinear coordinates equation (27) becomes 



(«,£) = £ L«, ^ (v.vm: (y,vj (28) 



'l L 'm = 1 ilL j„, 



An orthonormal basis for G(v p V 2 , K ,V m ,n) is given by, 
{djc^ 1 aL AdxJ" \ i l ,i 2 ,K i m e (1,2,K ,n)] and, 

(dx V] aL AdxJ™, cU,' aL Adx^) = 5 ; ,. L (5. . . (29) 

By extension, an orthonormal basis for G(v l5 K V V V 2 ,K ,V 2 ,K ,V m ,K v m ,n), where there 
are p ( . differentials of order V ( , and the v ; are distinct for i e {1, 2, K , m} , is given by, 

\dx v ) aL acU v ! A(k v ' aL Adjc v ' aL Adx v "; aL acU v '" I K , z' g (1,2,K 
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The coordinate indices assume all possible combinations with the constraint that 

1<^<!/ 2 <L < j pi <n, 1< 2 ^< 2 i 2 <L < 2 i pi <n,..., 1< J x < J 2 <L < J Pm <n. 



Let ae G(v p K ,v m ,n) = Yl Tl ^^ and P & G C"i> K >V L i> n ) = Tl W G ^ n 

7=1 1, M 

Then the exterior product of a and /5 obeys 



7=1 



) 

i=l 



aAjS = (-l)" jSa«. (30) 



The Hodge dual on G(v 1 ,V 1 , K ,V p n), in Cartesian coordinates, is constructed in the 
same way as in ordinary exterior algebra (see e.g., page 287 of [14] and page 108 of [21]). 
Specify a basis with an orientation {dx 1 v \ dx 2 ', K , doc v n l } for G(V[,n). The dxj' form an 
orthonormal basis on G(V[,n), (dxj 1 , dx^ ) = 8 ij . * denotes the Hodge dual, 

*(dx Vl aL Adx.') = dx v i ' aL Adx]' , (31) 

where (i,, K , i p ,j v K ,j„_ p ) is an even permutation of (1, K , n). In curvilinear 

n 

coordinates the Hodge dual of a = L , dy*' aL AdyJ 1 is also constructed just as it is 

i,L i p =l 

in exterior calculus but respecting the fractional coordinate transformation rules (see 
equation (57) of FDF I). Let J(v) = det(J/(x,y,v)) and 1/J(v) = det(J 7 ! (y,x,v)). Let 
£ J,h Jn denote the Levi Civita permutation symbol, i.e. £ J,h J "£ j L . = n\. Then, in index 



notation, 
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1 ^ £ hhi » 

(n-p)! jiL %Kv) 



~a = (-\) p{n - p) a 



If desired, ^jg(v) = ^jdet(g ij (x,y,v)) may be used in place of J(v) = det(J }(x,y, v)) . 



A dual for G(V P K V 15 V 2 ,K ,v 2 ,K ,V m ,K V m ,n) can be constructed by noting that, 



( p ^ 
G(v„K v„v 2 ,K ,v 2 ,K ,v m ,K v m ,n) = n ri G Cv f ,n) 

7=1 ^ i=l 



Then, 



*G(v 1 ,K ,v„v 2 ,K ,v 2 ,K ,v m ,K v m ,n) = n 



riG(v,,n) 

V 1=1 



*G(v 1 ,K ,v„v 2 ,K ,v 2 ,K ,v m ,K ,v m ,n) = n ri G ( v P n 



7=1 



v ,=1 y 



For example consider, 



a 



- E 



ft /? 



L ,• t ,• t ,• t ; ^ v , ! aL a dy v ) aL aJv"'" aL a dy v "; 

i'i L i' n L mh L „'„,„ - i<i y y m ! P , 
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In index notation, 



«= £ : i n 



1 £' 



i h L 1 < P1 L m ! i L m ! p„, . 



(38) 



a 



n(-ir (n ~ n) « 



(39) 



For the discussion of the inner product and Hodge dual for F(v,m,n) the restriction that 
V > is made. The inner product and Hodge dual for F(v,l,n) would be the same as for 
G(v,n). For F(v,2,n) , the basis elements are made up of two coordinate differentials. 



{dxf 1 Adxf | ij e {1,K ,n}y^ 2 > 0, 3 fl x + {l 2 = v] 



(40) 



Consider two arbitrary elements of F(v,2,n) 



n | p 

« = L J V ( a ij(yi> v - v i~) dx i i Adx T v< ) dv i> 



(41) 



P= £j„ V (A/V p V-V f X' A^P) dV t 



(42) 



The inner product on F(v,2,n) is, 



(a,jS) = £ J v (a .(v I .,v-v,.)A,-(v { ,v-v l .)) dv r 



1.7= 



(43) 
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In curvilinear coordinates equation (43) becomes, 



W)= £ ^(a^Cv^v-v^^Cv^v-v^Cy.v^Cy.v-v,) dv r (44) 



An inner product on F(v,m,n) would involve m-1 integrals and m summations. 

The problem with doing a Hodge dual on F(v,2,n) is that there is a subspace where the 

orders of the coordinate differentials are the same. That is G(v/2,v/2,n) . On this subspace 
the Hodge dual mapping will map into objects with (n - 2) coordinate differentials. On the 
rest of F(v,2,n) the Hodge dual mapping will map into objects with (2n - 2) coordinate 
differentials. To avoid this problem consider an arbitrary element of F(v,2,n) that has no 
component in G(v/2,v/2,n) . 



a = lim £ J l(a ij (v i ,v - vjdy? a dy v r v > ) dv t 



e->- 1,7=1 



Then, in index notation, the Hodge dual of a is, 



(45) 



a = lim 



V 

£->— 

2 



+ lim 

V fc 
£->— 



E 

n 

E 

'l,7l=l 



\2 



(n-1)! 



y'l L '„ C 7l L 7„ 



J K) J ( v - v J 



1 



\2 



(n-1)! 



g*l L '„ g7iL 7„ 



J(vJJ(v-vJ 



(46) 
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(47) 



The Hodge dual for F(v,m,n) is constructed in the same way, with the restriction of 

removing subspaces where the orders of the coordinate differentials are the same. Note that 
m < n is also required. 

3. Differential Properties. 

To begin this section, it is desired to determine the values of V and a for which the 
fractional Poincare identity is always true. V being the order of the exterior fractional 
differintegral and a being the object that is acted upon. 



d v d v (a) = 



(48) 



Let ae G(v l5 K ,Vj,V 2 ,K ,V 2 ,K ,V m ,K ,V m ,n) and denote the elements of an orthonormal 

N 

basis by a, . To be specific, a = ^ccp, , N = 



i=i 



n 



n 

t 



n 

\P2J 



\FmJ 



<7,& |ck|'! aL Adx^l aL Adx v ™ aL Adx^ I x i xn i 2 K m i p e (1,2,K and thepjL p n 

are defined in equations (22) and (23). Substituting this into equation (48) and expanding 
the fractional exterior differintegral gives, 



N n 



dM , '(a) = Id l 'I(k;A fl; D;(a / a / ). 

i=i j=i 



(49) 



The product rule for fractional differintegrals is used to evaluate D v x (cc^,) (equation (13) 
of section 1). 
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The second factor on the right-hand side will be zero for all k values except k = 0. Thus, 

.pliap^pia,)*,, (51) 

and equation (49) becomes, 

N n 

d"d"(a) = Ed"E( B D;a,)dx;A<T I . (52) 

1=1 7=1 

The second fractional exterior differintegral is expanded and reduced in the same manner 
yielding 

d v d>) = LL( fl Dl ( ( fl D: a,))(dx: Ad*; act,). (53) 

i=i j,k=i 

The factor dx v k Adjc^ act, is antisymmetric in the k and j indices. The fractional 

differintegral operators are linear and with respect to different coordinates and are thus 
symmetric on the k and j indices. Hence, each term in the above summation is zero. Thus 
the fractional Poincare lemma is true for all values of V e C and for all 
ae G(v p K v v v 2 ,K ,v 2 ,K ,v m ,K v m ,n) (assuming that the components of a are 
differintegrable to the appropriate orders). Recall that, F(v,m,n) is the union of all the 
G(v 1; v 2 , K ,v m ,n) i.e. 
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F(v,m,n)= U G(v 1 ,v 2 , K ,v m ,n). (54) 



Hence, 



d v d v (a) = VveCandVae F(v,m,n) . (55) 

In FDF I the coordinate transformation rules were worked out for fractional differential 
forms. For fractional integral forms the coordinate transformation rule is less compact. 
The reason for this is that fractional derivatives (Re(v) > 0) have a non-trivial kernel (see 
equation (31) of FDF I) while fractional integrals (Re(v) < 0) do not. To begin, assume 
that the Cartesian coordinates, {*,.}, can be written smoothly in terms of the curvilinear 

coordinates, {y ( . } , and note the expression for the exterior fractional differintegral in the two 
coordinate systems. 

x i =x i (y) (56) 

d v = X>r a D; (57) 

i-l 

d v = I>L I)!: (58) 

7=1 

The goal is to find a coordinate transformation matrix, J ! J (x,y,v), that will express the 
dxj in terms of the dy, v . 
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7=1 



(59) 



To construct the matrix for the integral order forms let k e {1,K ,n}. Then, the following 
will generate n linear equations for the dxj . 

I>r a , % m = i>; ~ aj d; (60) 

i=l j=l 

Cramer's rule can now be used to find the matrix for the coordinate transformation. Define 
the n x n matrix A to be, 

A,* = a p\{x k ), (61) 

and a 1 x n column vector, 

b* = Ldv^,, D: (.v i( v)). (62) 

7=1 

Equation (60) can now be written as, 

X>, V A ( .,= b k . (63) 

The cbc, v can now be solved for in terms of the dy v .. Define , A to be the matrix A with 
column number i replaced with the column vector b^. Then, via Cramer's rule, 
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(64) 



Equation (64) defines the coordinate transformation matrix for fractional integral forms. 



< = L^C-W) (65) 



7=1 



As an example, consider the transformation from Cartesian to polar coordinates with the 
order v = -1 and the initial point for the fractional integrals is the origin. 



2tan(0)-l , 2(tan 2 (0) + 2) 

dx" 1 = W dr-' - — '-dO- 1 (66) 

3 3rsin(0) 



2cotf0)-l , 2(cot 2 (0) + 2) 
3 3rcos(0) 



To construct a covariant fractional derivative the restriction Re(v) > is made. A 
covariant fractional derivative must have the following properties. 

lim ;V„ =V, (68) 



V->1 



.Km 'V^DJ (69) 

{y}^{x} 



£./;(v.-.v);v v =;v. (7(» 
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Equation (68) is merely a statement that we get the covariant derivative that we are familiar 
with from differential geometry; i.e. covariantly constant metric and parallel transport. 
Please note that the metric need only be covariantly constant for the case of derivatives of 
order 1 and for a background space of order 1 . This is because parallel transport is tied to 
the notion of the tangent vector. Equation (69) says that when the curvilinear coordinates 
become Cartesian the covariant fractional derivative becomes an ordinary fractional 
derivative. Just as in differential geometry the covariant derivative becomes a partial 
derivative when space becomes flat and the coordinates Cartesian. Equation (70) is just the 
statement that the fractional covariant derivative transform as a vector. 

From FDF I the following transformation rules are available, 

aPl^'jiy^nPlr (71) 
i=i 

V*(*) = XX(y,x,v)V a (y). (72) 
The covariant derivative of a vector must transform as, 

Ij',(v,x,v)j";(v,^,v)(, ; ;v yi v m (v))= p:y k ( X ). (73) 

i,m=l 

Now substitute the known transformation rules into equation (73). 

tj'^v)!'^ (74) 

i,m=l a,i-\ 
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n n 

Multiply equation (74) by ^J7(x,;y,v) and recall that, j(y,x,v)i h J (x,y,v) = d' h gives, 



7=1 



£n(y,x,v)( a ;V y V m (y)) = £ 56 D;(j a ,();,x,v)V a (3;)). (75) 



m=l «-l 



Multiply equation (75) by ^J, k (x,y,v) , simplify again, and the fractional covariant 



derivative is seen to be 



lV v V / ();)=i;j / i (x,v,v)^D; 6 (j;( J ,^,v)V fl (v)). (76) 



a* 

*,a=l 



If V — > 1 then the expression reduces to, 



k,a=\ k,a=\ 



Equation (77) is the usual expression for a covariant derivative, written in terms of the 
coordinate transformation matrix. 

Now expand equation (76) using the product rule for fractional derivatives 

fi :^V,(y) = £ J^^.v^Id^^J;^^)). (78) 



k,a=\ 5=0 



If the first term of the infinite series is separated out the connection can be isolated. 
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5 rV y V,(y)=£V(*,y,v)|( 5 ^^^ 

(79) 



k,a=l 



^^00=^,00+ £ h k (w)t[ V X Drv.(y)Rj- 4 (y,*,v)) (80) 

*,a=l .s=l V / 



This is somewhat more complicated than the usual connection coefficients that arise when 
constructing a covariant derivative. So, adopt the following notation and refer to the 
connection as a connection functional (due to the integrals involved in the sum). 



,IV y V,(v)= a -D;v i( y) +5 :7,(V) 



(81) 



k,a= 



5=] 



*>»(V) = £ J,*(w)£ V kj\(v,x,v)( 5) D;;x(y))) 



(82) 



As a cautionary note, recall that these objects are defined on a Euclidean space, i.e. there 
is no curvature. For a curved manifold there would be difficult issues concerning the non- 
local nature of fractional derivatives. Another difficult problem is to find an expression for 
a covariant derivative that is a different order than the background space. How does one 

construct £V y on a space where the vectors transform under J%(y,z,v)7 The obvious path 

to approach this problem would be to use the chain rule. However, the chain rule in 
fractional calculus is computationally difficult to use for an arbitrary order (see pages 97 
and 98 of [21]). For an alternate discussion of the fractional covariant derivative see 
Kobelev [11]. 
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IV. Matrix Order Forms. 



In the appendix fractional derivatives are extended to matrix order derivatives. Matrix 
order derivatives are shown to be well defined for all square matrices over the complex 

numbers ( c mxm ). Consider the definition of a fractional exterior differintegral. 



(83) 



To construct a matrix order form, replace the parameter ve C by a matrix A e C' nxm . This 
can be expressed directly or, using the spectral theorem (see page 517 of [15]). 



d A = £d* A a D A =£d* A P 



D 



o 



p-i 



(84) 



(85) 



Where A = PDP D is a diagonal matrix whose entries are the eigenvalues of A, 

k 

A = Eg k is the number of unrepeated eigenvalues, A,, are the eigenvalues, and the G, 

!=1 

are the spectral projectors. 



To develop an understanding of these new objects, consider the coordinate 
transformation rules for matrix order forms. 
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(86) 



For now assume that A is diagonalizable, the Jordan case will be dealt with in a later paper. 
Express the Cartesian coordinates in terms of the curvilinear coordinates, x i = x^y), and 
compute the exterior fractional differintegral in the two coordinate systems. 



d A (x k ) = d A (x k (y)) 



(87) 



Id*, A a D^)=I>^ D A (x k (y)) 

i=i j=\ 



(88) 



To find the coordinate transformation matrix equation (88) needs to be solved for the dxf. 
a D A (x k ) and s D A (x k (y)) form nxn matrices as i,j,k& {1, K ,n} . Hence equation 

(88) represents n 2 equations for the n dx A . This is an over determined system, hence, dx A 

cannot be viewed as a single object. The simplest thing to do is to view dxj as an object 
continually depending on the parameter v, and then write (see page 526 of [15]), 



(dx^ 



dx A =P 







dx> 



P 1 . 



(89) 



With this assumption consider equation (88) 
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n 



dx^ 







dbc> 



7=1 



o 



P-'sDj^W). (90) 



Equation (90) can be simplified if a and a D^ 1 are expressed in terms of their diagonal 
matrices and all factors of P" 1 and P are either multiplied together or factored out. 



i=l 



dx^ D^(x t ) 







"•~! a i X; k' J 



=E 

7=1 



dyf a p^(x k (y)) 







dy^ 5 p^(x k (y)) 



(91) 



This can be expressed as a single equation, 



(92) 



7=1 



where /e {1,K ,m}. A coordinate transformation matrix can be constructed using 
Cramer's rule as was done in section 4. If all the eigenvalues are such that Re(A / ) > then 
the construction from FDF I can be used. In either case a prescription is available for the 
construction of the coordinate transformation matrix for matrix order forms. 



dx 







dx 



=E 

7=1 



dyf.r^y) 







(93) 



If P and P 1 are inserted in the appropriate places equation (93) becomes, 
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(94) 



where, 



rj(A,x,y) = P 







P 1 . 



(95) 



With the coordinate transformation matrix found a metric and line element can be 
constructed in the same manor as in FDF I. A matrix order covariant derivative can also be 
constructed provided that attention is restricted to matrix orders that are positive definite. 



aXViO'HDjv.CyH^y.CV) 



(96) 



y»00 = £ V(*,y,A)£ P D^V a (v))«J\(j,x,A)) 

k,a=l s=l V ) 



(97) 



J, k (A,x,y) is defined in equation (95), the other components of equation (97) are defined in 



equations (98), (99), and (100). 



"V 

V s J 







K S J) 



(98) 
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D A_lI = P 

i y'b 



yb 







P" 1 



(99) 



d sl = 



yb 







<9 S 



(100) 



I is the m x m identity matrix. Equations (95)-(99) can also be expressed using the spectral 
theorem. 



The Poincare lemma also holds true for matrix order exterior differintegrals. Let 
ae G(V[,K v v v 2 ,K ,v 2 ,K ,v k ,K v k ,n) and allow the coefficients of a to possibly be 

matrix valued, then d A d A a can be written as, 



d A d A a = £ £p 

,=i j=i 







( JCij Xj CI j Xj j 



P"'a. (101) 



Each component of the right-hand side of equation (101) is now equivalent to the problem 
considered in section 3. Hence, d A d A a = for any A e (j mxm that is diagonalizable and 
for any a e F(v,m,n) (assuming that the components of a are differintegrable to the 
appropriate orders). For matrices that are only Jordan diagonalizable the result does not 
hold. This is because the upper triangular Jordan blocks do not commute. 
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V. Conclusion. 



In this paper the results of FDF I were extended to include an inner product, Hodge 
dual, and covariant derivative for fractional differential forms. The connection for covariant 
fractional order derivatives was found to be somewhat more complex than for covariant 
derivatives of order one. The resulting formula does reduce to the usual formula from 
differential geometry. The notion of a fractional differential form was also extended to 
fractional integral and matrix order forms. Coordinate transformation rules were worked 
out for both of these objects with a specific example presented for integral order forms. 
Matrix order forms were constructed by combining matrix order fractional calculus (see the 
appendix) with exterior derivatives. This was done in the same way that fractional order 
derivatives were combined with exterior derivatives in FDF I. The Poincare lemma was 
found to be true for all orders of exterior differintegrals provided that the order of the 
differintegral was a diagonalizable matrix or any complex number. 

One final note, the Riemann-Liouville fractional differintegral was used for the 
construction of fractional forms. There is no reason that the Griinwald-Letnikov, Weyl, or 
the Caputo definition of fractional differintegrals could not be used to define fractional 
forms. See references [17] and [21] for the definitions of these other fractional 
differintegrals. 

Acknowledgment: The authors would like to thank P. Dorcey for helpful comments and a 
critical reading of the paper and Jean Krisch for asking the right questions to get all of this 
started. 
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Appendix: Matrix Order Differintegration 

In this appendix the formulae from the Riemann-Liouville formulation of fractional 

calculus are adapted to matrix order for any matrix A e C nxn (nxn square matrices over 
the complex numbers). A brief review of matrix properties is provided to set notation and to 
provide a convenient reference for the reader. 

The concepts of integration and differentiation have been expanded many times 
throughout the development of calculus. Almost immediately after the formulation of 
classical calculus by Leibnitz and Newton the question of half order derivatives arose (see 
[16] and [17] for historical reviews). Derivatives of complex and purely imaginary order 
were later developed (see e.g. [13] and references there in). Phillips considered fractional 
order derivatives of matrix-valued functions with respect to their matrix arguments, 
fractional matrix calculus. These operations were found to be well defined and 
subsequently applied to econometric distribution theory (see [18]-[20]). 

To develop the notion of matrix order derivatives and integrals (differintegrals) it is 
useful to review the properties of functions with matrix order arguments. Three different 
means of computing functions with matrix order arguments will be used in this paper. The 
method chosen will be determined by the properties of the matrix being used and the 
specific application being considered. 

A given non-zero matrix A e c mxm can be placed into one or more of the following 
sets; Jordan block diagonalizable, diagonalizable, and/or normal. These sets will be denoted 
by B m , D m , and N m , respectively. Note that N m c D m c B m and every non-zero matrix 

A e C mxm is at least Jordan block diagonalizable. Recall that the set of normal matrices 
contains the following subsets: real symmetric, real skew-symmetric, positive and negative 



27 



definite, Hermitian, skew-Hermitian, orthogonal, and unitary (see page 548 of [15] for 
further details). If a matrix A is normal then there exists a unitary matrix U such that 

A = UDU*. (102) 

Where D is a diagonal matrix with the eigenvalues of A as the entries and * denotes 
conjugate transpose. If A is real and normal then equation (102) reduces to 

A = ODO T (103) 

and O is an orthogonal matrix. If a matrix A is diagonalizable then there exists an invertible 
matrix P such that 

A = PDP-'. (104) 

Alternatively, if A is diagonalizable the spectral theorem (see page 517 of [15]) can be used 
to find another representation for A. If A e D m then there exists a set of matrices 
{G p K , Gj such that, 

A = £G,A,, (105) 

i=i 

where {A,, K , X k } are the unrepeated eigenvalues of A. The matrices {G 15 K , G,.} have 
the following properties; 

G,G. =0 for i* j, (106) 
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G,G, =G, , 



(107) 



I = £g, . (108) 
If a matrix is Jordan block diagonalizable, then there is an invertible matrix P such that 



A = PJP" 1 = P 







•I 



o KK) 



(109) 



Where \, K , ^ are the distinct eigenvalues of A ( < m ), and J(Aj) is the Jordan 
segment for the eigenvalue "k- . 



[\ 1 





1 

X; 



(110) 
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There are several ways to define a matrix function. Let A e N m with eigenvalues 
\, K , A m (they need not be distinct) and let g(z) be a function that is defined for all the 
eigenvalues of A. Then g(A) can be expressed as 



to 



g(A) = U 







u*. 



(in) 



If A is real replace U by O and U by O . If A is diagonalizable then equation (111) 
would be written as, 
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g(A) = P 







gaj 



rn ' J 



(112) 



The spectral representation can also be used (in either case). 

g(A)=LG,.ga,.) 

Where \, K , \ are the distinct eigenvalues. For matrices that are Jordan block 
diagonalizable but not in D m the situation is somewhat more complicated. 



(113) 



(g(W) 



g(A) = P 







g(W), 



(114) 



Where, 



g(W) = 



f g(A,) g'a ; ) g"(A,.)/2! L gC-'^O/a-D! 
g(A,) g'(A,.) L g('- 2 )(A,.)/(/-2)! 

M 



g(^) 



g(^) 



(115) 



g(J(A,)) is an / x / upper triangular matrix. Note that for this case the function need be at 
least differentiable of order / - 1 at X i . For further details see [ 1 5] . 
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Using equations (4) and (5) of section 2 an analytic function can be defined (see page 
49 of [17]) for a given f(x) that is suitably differentiable. 



S(A)= D A A f(*) 



1: 



Re (A) < 0, 



1'- 



Re (A) > 



n > Re(A) (n is whole) 



(116) 



Let A g N m then the differintegral of order A is given by, 



D 1 ' 

a x 



X>t = U 



o 



a x J 



LP 



(117) 



If A is real and normal then, 



X>Z =0 







D 



o 1 



(118) 



If A e D m then 



f J N 

D 1 

a x 







D 

a x J 



p-i 



(119) 
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Any of the above three cases, equations (1 17), (118), or (1 19), can be expressed using the 
spectral theorem. 



(120) 



Where the G, are the matrices as given by the spectral theorem and the sum is over the 
unrepeated eigenvalues. 

For matrices that are Jordan block diagonalizable but not diagonalizable some notation 
must be introduced. Recall equation (116). 



g a)= a D x x f( X ) =• 



-LJ'J^. Re (A)<0, 

r(-A) K (x-^ + " y ' 



dx" 



X-n+l 



Y{n-X) Ja {x-%) 



Re (A) > 



n > Re(A) (n is whole) 



(121) 



Derivatives of g(X) with respect to X will appear in the expression for the Jordan matrices. 



iLlJ'JS^, Re(A)<0, 



d k d n 



^1' 

1 \ *> n 



X-n+l 



Y{n-X) Ja {x-%) 



Re(A) > 



n > Re(A) (n is whole) 



(122) 



In equation (122) k is restricted to be a whole number. The J) x can be expressed in terms 
of a D^. For example, if k = 1 the following is obtained. 
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d_ 



ln(-A) + ln(*) [D, A /(x)-£ aD " (x Y (X)) Re(A)<0 (123) 



i D V(jc) = -— 



7=1 



Re (A) > 



n > Re(A) (n is whole) 



(124) 



Similar expressions are available for the higher order derivatives of a D x . With this notation 
matrix order differintegrals for A e B m are denoted by 



D A =P 

a x 



D 

a x 







D 



7U t ) 



P" 1 



(125) 



where, 



(°D X > *dV1! L 



D 



,M ) _ 



D A, 

a x 



L 




:Dj/(/-2)! 



M 



(126) 



To determine some of the properties of matrix order differintegrals consider the 
composition of two matrix order differintegrals. Let A and B e D m such that A = PDP" 1 

and B = QEQ" 1 where D and E are diagonal matrices with the eigenvalues of A and B as 
entries. Denote the eigenvalues of A as A, and the eigenvalues of B as p r Then, 
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( D \ 



D A D S =P 

a x a x 







P _1 Q 



a x 







D p < 

a x J 



Q" 1 . (127) 



To simplify this denote by R = P 'Q and the components of R as R, v . Equation (127) now 
has the form, 



D A D B =P [R D p 1 Q" 1 . 

a x a x y ij a x a 1 J ^ 



(128) 



There is no sum over the repeated indices, they merely denote the components of the matrix 

between P and Q" 1 . Alternatively the spectral theorem can be used to obtain another 
representation for the right hand side of equation (127). 



k I 



d a d b =Vg t> k Tn D P =V Y"gh d Pj 

a x a x • a x L~ • j a x L~ 1^ i j a x a x 



K t 

Where A = £g, X i and B = ^H , p. 



Equation (4) carries over to the matrix order case. 



(129) 



~\m ~\m k 

TvD A f(x) = — LG lB Djf(x) 

ox ox (=1 



(130) 



k ~\m 

r£ dx ma x 



(131) 
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(132) 



(=i 



= Df m f(^) 



(133) 



I is the identity matrix, m is a whole number, and f(x) may be a scalar, vector, or matrix 
order function. 

Equation (7) carries over in the following sense, Let A be a matrix such that Re(A,) > 0. 
Now consider equation (127) with B replaced by -A. 



D A D A =P 

a x a x 



o 



D" 







P" 1 (134) 



a x a x 



D A D A =P 



a x a x 



o 



a x a x 



P" 1 



(135) 



Which, by equation (7), reduces to the identity operator. 

Suppose that A and B commute (and are diagonalizable), then A and B can be 
diagonalized by the same matrix, say P. Equation (127) is now, 



a x a x 



D D*=P 

a x a x 



o 



a x a x 



P 1 . 



(136) 
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A further simplification of equation (136) will occur using equations (5), (7), (8), (9), or 
(10) as the signs of the eigenvalues dictate. For example, if the eigenvalues of A and B are 
such that Re(A ; ) < and Re(p ; ) < Othen equation (136) becomes, 



D A D B =P 

a x a x 







T\Kn+Pm 

ax J 



-1 _ p»A+B 

ax 



(137) 



If the eigenvalues of matrices A and B are such that Re(A, ) < and Re(p,) < but not 
commuting then equation (127) is, 



D A D B =P [r D A +Pj 1 Q" 

a x a x [ lj a x J ^ 



(138) 



Or, using the spectral theorem, 



k I 



a D, A fl D!=IlG,i/ Jfl D 



(139) 



Let A, B g N m with A = PDP T and B = QEQ T where D and E are diagonal. 
Denote the eigenvalues of A by X i and the eigenvalues of B by p r Now consider the 
transpose of equation (127). 



( D A D B ) T = Q 



a ^ x 









a x J 



Q T P 



D 

a x 



o 







D , 

a x j 



(140) 
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= D B D A 

a x a x 



(141) 



As a final result consider the determinant of equation (119). 



det( fl D A ) = det(P) 



V i=l ) 



i=l 



The right hand side of equation (142) is a sequential fractional derivative (see e.g. [16] and 
[17]). Thus a sequential factional derivative may be viewed as the determinant of a matrix 
order fractional derivative. Additionally if the eigenvalues of A are such that Re(A, ) < 
then equation (142) can be simplified to give, 

det(X) = fl D A Tr(A) (143) 
where Tr(A) is the trace of the matrix A, i.e. the sum of the eigenvalues. 
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